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FIG. 1. Optical image and sketch of 10 x 10 square array. 

The junction area is 9 fim^ and their critical cur- 
rent density is about l.lKA/cm^ for square array, and 
O.lKA/cm^ for triangular array, with a typical spread 
of junction parameters of about 5%. The square ar- 
rays consist of 10 columns by 10 rows and have the cell 
area S = 38 ^m^. The triangular arrays are made of 12 
columns by 12 rows and have S — 160 /zm^. All arrays 



FIG. 3. (a) Enlargement of digitally stored I-V curves in 
the region between zero and 4Vg; (b-e) LTSLM images of 
resistive states of the array at the bias points close to the 

By using the LTSLM, we have systematically imaged 
different resistive configurations biased at various fine re- 
sistive branches of the I-V curves. The experimental pro- 
cedure is the following: the number n of switched rows 
is fixed by biasing the array at voltages V ~ nVg. Then 
LTSLM images of the sample are recorded at a constant 
bias value. Typical images from the square array are 
shown in Fig. 3(b)-(e). The black spots correspond to 
the junctions that are in the resistive state. The junc- 
tions that are in the superconducting state do not appear 
on the array image. As expected, the images show that 
the number of rows which are in the resistive state is 
equal to the number n of gap voltages selected by the 
bias point. To trap a different resistive configuration, we 
always increased the bias current / above I^"^^ , the ar- 
ray critical current, and then reduced it to the low level. 

Similarly to the experimental results of Ref. H, we 



shown in Fig. 5. Resistive states with both single and 
double meanders have been trapped (Fig. 5). The ap- 
pearance of the meandering is qualitatively similar to the 
experimental images shown in Fig. 3(b-e). The simula- 
tions have been carried out in the presence of finite mag- 
netic field in order to prevent the simultaneous switching 
of all the array rows It leads to the decrease of the 
critical current and the appearance of the flux-flow region 
where I-V curves are highly non linear (Fig. 5). 

We have found that, in the absence of disorder, the 
probability of the appearance of a meander in two row 
array is about 1.8% per horizontal junction. Surprisingly, 
in the presence of specially introduced disorder up to 20% 
the numerical simulations show no increase of the rate of 
appearance and no preferred position of the meanders. 
Due to this fact, we conclude that the broken symmetry 
row switching appears due to an intrinsic instability 



FIG. 4. LTSLM images of more complicated dynamic 
states observed at the bias points close to 7Vg voltage region. 

In the following, we present numerical simulations of 
the dynamics of underdamped arrays. Calculations were 
performed using the resistively shunted model for Joseph- 
son junctions and the usual analysis for superconductive 
loops with only self-inductances taken into account, i.e. 
the Nakajima-Sawada equations For a discussion of 
these equations and the application limits see Refs. 1^,1). 
The phase configuration of a row switched state was used 
as initial condition for the next run. 

We simulated large (10 x 10) and small (2 x 10) square 
arrays with various parameter (3^ between 0.5 and 2. In 
all cases, the simulations well reproduce the branching of 
I-V curves and meandering character of the finite voltage 
paths. Moreover, in order to check the influence of self- 
inductance effects the simulation were carried out for the 
arrays with a small parameter Pl = 0.1. A simulated 
I-V curve with three different branches and the corre- 
sponding voltage distributions for the 2 x 10 array are 
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FIG. 5. Numerically simulated I-V characteristics of a two 
row array and the corresponding voltage paths at 7 = 0.4. Jc 
is the single junction critical current. 

of the superconducting state, i.e. a change of the ini- 
tial conditions can lead to the appearance of a different 
percolative voltage path. The numerical analysis of the 
two row array also allows to find out some peculiarities of 
the new dynamic states. When the junctions switched to 
the resistive state form a straight line in the top or bot- 
tom row (no meandering), the array has the highest resis- 
tance, Rn/ {N-\-l) (triangle symbols in Fig. 5), where TV is 
the total number of cells in the row and i?„ is the normal 
resistance of one junction. In the case when the volt- 
age path with one meandering occurs, the array shows a 
lower resistance, i?„/ {N -\- 2) (square. Fig. 5). Now there 
is one more junction connected in parallel to the resistive 
path, with respect to the case of no meandering. For the 
same reason, the meanders across two horizontal junc- 
tions corresponds, to the resistance i?„/(iV -t- 3) (circle. 
Fig. 5). The presence of meanders shifts the array I-V 
curve to the left and this resistance scaling explains the 
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experimentally observed fine branching of the I-V curves. 

We have found that the above dynamic states are sta- 
ble in a wide range of voltages, i.e. we have not observed 
direct switching between branches. This is due to the 
too large energy necessary for the junction capacitance 
recharge in underdamped arrays [Q. So, only switch- 
ing from the superconducting state allows to trap these 
dynamic states. 

Our analysis using the KirchhofF's current laws shows 
that, in the case of meandering, the dc superconducting 
current flows via the horizontal junctions. Moreover, this 
current increases from the row boundary to the cell where 
the meandering occurs. This current configuration is dif- 
ferent from that without meanders. In the latter case the 
horizontal junctions are not active and dc superconduct- 
ing current flows only via the vertical junctions. 
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FIG. 6. Sketches of metastable superconducting and re- 
sistive states for a small array with two plaquettes containing 
three junctions per cell: (a) mesh currents are equal to zero, 
(b) mesh currents are different from zero. Thick solid lines 
show the parts of the array switched to resistive states, and 
numbers are the local dc currents normalized to Ic. 

All these features can be explained, at least quali- 
tatively, by the presence of metastable superconduct- 
ing states in 2D Josephson arrays. We analyze these 
states for a simple case of two plaquettes with 5 junc- 
tions (Fig. 6). In the absence of magnetic field, the most 
stable superconducting state is the one with all mesh 
currents equal zero. At the bias current / = 2/c the 
vertical junctions switch from superconducting state to 
resistive state with straight resistive rows (dc mesh cur- 
rents are still zero) (Fig. 6(a)). However, this system has 
another peculiar state with finite values of mesh currents, 
as shown in Fig. 6(b). At I = Ic this system can support 
a superconducting state with the Josephson phases of the 
horizontal junction and one vertical junction per row be- 
ing equal to 7r/2. The Josephson phases of the remaining 
vertical junctions are equal to tt. This metastable state 
supports finite mesh currents /c/2 and switches to the dy- 



namic (resistive) state with a meander via the horizontal 
junction. The mesh currents decrease to the value of /c/6 
and change sign during the switching process, which re- 
sults in a stable dynamic state with meandering. Thus, 
the dynamic states with meandered dc voltage drop are 
" fingerprints" of metastable superconducting states with 
non zero mesh currents and active horizontal junctions. 
The state appears only when there are at least three junc- 
tions per cell and it exists also in the limit of zero lin- 
ear self-inductance and zero mag netic field @. That 
is different from the case of single- and double-junction 
SQUIDs, where metastable superconducting states ap- 
pear only in the limit of large self- inductance |l2| . Since 
the above discussed metastable state exists in the limit of 
zero inductance, the small arrays with few plaquettes can 
be suitable system to observe the effect of macroscopic 
quantum coherence at low temperatures |r^ , [l5|| . 

It appears that also large 2D arrays can provide var- 
ious metastable superconducting states that lead to the 
broken symmetry in row switching process. Theoretical 
and experimental investigation of these metastable super- 
conducting states and their dependence on the externally 
applied magnetic field are in progress. 
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